The Floquet Method for PT-symmetric Periodic Potentials 
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By the general theory of PT-symmetric quantum systems, their energy levels 
are either real or occur in complex-conjugate pairs, which implies that the secular 
equation must be real. However, for periodic potentials it is by no means clear that 
the secular equation arising in the Floquet method is indeed real, since it involves 
two linearly independent solutions of the Schrodinger equation. In this brief note we 
elucidate how that reality can be established. 
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The study of quantum mechanical Hamiltonians that are PT-symmetric but not 
Hermitian[l|-j6[ has recently found an unexpected application in classical optics0j-[l2j, due 
to the fact that in the paraxial approximation the equation of propagation of an electromag- 
netic wave in a medium is formally identical to the Schrodinger equation, but with different 
interpretations for the symbols appearing therein. The equation of propagation takes the 



o f form 

A '' % = -(£ + v<.*)H, (i) 



dz \dx 

where ip(x,z) represents the envelope function of the amplitude of the electric field, z is a 
scaled propagation distance, and V(x) is the optical potential, proportional to the variation 
in the refractive index of the material through which the wave is passing. That is, V(x) cx 
v(x), where n = n (l + v(x)), with \v\ <C 1. A complex v corresponds to a complex refractive 
■ index, whose imaginary part represents either loss or gain. In principle the loss and gain 

regions can be carefully configured so that v is PT symmetric, that is v*(x) = v(—x). 

In the optical context one is often interested in periodic PT-symmetric potentials, which 
have a number of very interesting properties [3]- 12]. Such a potential V(x), whose period 



we can take as ir, without loss of generality, satisfies the two conditions V*(— x) = V(x) 
V(x + 7r). For a periodic potential we are interested in finding the Bloch solutions, which 
are solutions of the time-independent Schrodinger equation 



X. 
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+ V(x))Mx) = EMx) ( 2 ) 



dx 

with the periodicity property if>k{% + tt) = e lk7T ipk{x)- 

The Floquet method consists of expressing ipk(x) in terms of two solutions, u\(x) and 
U2(x), of Eq. (|2}, with initial conditions 

ui(0) = 1, ui(0) = , 

u 2 (0) = 0, u' 2 (0) = 1 . (3) 
Then ipk{%) is written as the superposition 

^k{x) = c k ui(x) + d k u 2 {x). (4) 
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Imposing the conditions ipki 7 ^) = e tk7T ilj(0) and ip' k {^) = e ij)'(0) and exploiting the invari- 
ance of the Wronskian W{ui,u 2 ) one arrives at the secular equation 

cos kn = A = - (ui(n) + u' 2 (n)) . (5) 

In the Hermitian situation both Ui{tt) and u 2 {tt) are real, and the equation for k gives real 
solutions (bands) when |A| < 1. However, in the non- Hermitian, PT-symmetric, situation 
it is not at all obvious that A is real, since that implies a relation between u\{tt) and u' 2 (7r), 
even though u\(x) and U2(x) are linearly independent solutions of Eq. (T5]). It is that problem 
that we wish to address in the present note. In fact we will show that u 2 (tt) = ul(7r). 

The clue to relating ui(ir) and u 2 (ti) comes from considering a half-period shift, namely 
x = z + 7r/2. We write <p(z) = ip(z + 7r/2) and U(z) = V(z + tt/2). Then <p(z) satisfies the 
Schrodinger equation 



— + U{z)) Vk {z) = Ei, k {z). (6) 



a 2 

dz 

The crucial point is that because of the periodicity and PT-symmetry of V(x) the new 
potential U(z) is also PT-symmetric. Thus U(— z) = V(—z + 7r/2) = V(—z — n/2) = 
V*{z + ir/2) = U*(z). 

Now we can express the Floquet functions U\(x), U2(x) in terms of Floquet functions 
Vi(z), v 2 (z) of the transformed equation ([H]), satisfying 

v l (0) = l, v[(0) = 0, 

v 2 (0) = 0, ^(0) = 1 . (7) 
It is easily seen that the relation is 

Ul ( x ) = v 2 (-7r/2)v 1 (z)-v' 1 (-7r/2)v 2 (z) 

u 2 (x) = -V2(-ir/2)vi(z) +v 1 (-n/2)v 2 (z), (8) 

in order to satisfy the initial conditions on ui(x), U2{x). 
So 

Wi(tt) = v' 2 (-7i/2)v 1 (7i/2)-v' 1 (-n/2)v 2 (n/2) 
u[(n) = v 2 (-ir/2)v[(n/2) -v[(-it/2)v 2 (it/2) 
u 2 (ti) = -v 2 (-7T/2)v 1 (n/2)+v 1 {-7T/2)v2(7i/2) 

u' 2 (ir) = -«2(-7r/2)ui(7r/2)+i;i(-7r/2)^(7r/2) . (9) 
But, because of PT-symmetry, 

v 1 (-tt/2) = Mtt/2))* 
v[(-n/2) = -(v[(7r/2)r 
v 2 (-n/2) = -M*/2)y 

v' 2 (-rr/2) = (v' 2 (n/2)y . (10) 

Hence, indeed, ui(n) = {u 2 {n))* , and moreover, u[(rr) and u 2 (n) are real. The statement 
ui(tt) = {u' 2 (tt))* is in fact the PT-generalization of the relation Ui(it) = u' 2 (tc) implied 
without proof by Eq. (20.3.10) of Ref. [13j for the Hermitian case of the Mathieu equation, 
where V(x) = cos(2x). 
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If we wish, we may express everything in terms of u\, u 2 because from Eq. 



Ui(tt/2) 

ui(7T/2) 

u 2 (n/2) 
u 2 (tt/2) 



Vi(-7t/2) 



Hence 



Ul (n) = K(tt))* = Wl ( 7 r/2)K(7r/2))* + U ' 1 (7r/2)(« 2 ( 7 r/2)) s 



(11) 



(12) 



which is the PT-generalization of a relation implied by Eq. (20.3.11) of Ref. [13j after the 
use of the invariance of the Wronskian. 
Similarly 



u[(7t) = 2Re(w* (7r/2)ui(7r/2)) 
« 2 (vr) = 2Re(M;(7r/2)M' 2 (7r/2)) 



(13) 
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